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ABSTRACT 

The ultraviolet spectrum (1145-1720 A) of the distant quasar Q 0302-003 
[z = 3.286) was observed at 1.8 A resolution with the Space Telescope Imaging 
Spectrograph aboard the Hubble Space Telescope. A total integration time 
of 23,280 s was obtained. The spectrum clearly delineates the Gunn-Peterson 
He II absorption trough, produced by He II Lya along the line of sight over the 
redshift range z = 2.78 — 3.28. Its interpretation was facilitated by modeling 
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based on Keck HIRES spectra of the H I Lya forest (provided by A. Songaila 
and by M. Ranch and W. Sargent). We find that near the quasar, He II Lya 
absorption is produced by discrete clouds, with no significant diffuse gas; this 
is attributed to a He II "proximity effect" in which the quasar fully ionizes 
He in the diffuse intergalactic medium, but not the He in denser clouds. By 
two different methods we calculate that the average He II Lya opacity at 
2; 3.15 is r > 4.8. In the Dobrzycki-Bechtold void in the H I Lya forest 
near z — 3.18, the average He II opacity r = 4.47lo;33- Such large opacities 
require the presence of a diffuse gas component as well as a soft UV background 
spectrum, whose softness parameter, defined as the ratio of the photo-ionization 
rate in H I over the one in He II S" = ^m/^ucii — §00) indicating a significant 
stellar contribution. At z = 3.05, there is a distinct region of high He II Lya 
transmission which most likely arises in a region where helium is doubly ionized 
by a discrete local source, quite possibly an AGN. At redshifts z < 2.87, the 
He II Lya opacity detected by STIS, r = 1.88, is significantly lower than at 
z > 3. Such a reduction in opacity is consistent with Songaila's (1998) report 
that the hardness of the UV background spectrum increases rapidly from z — 3 
to z^ 2.9. 

Subject headings: cosmology: observations - galaxies: quasars: absorption lines 
- galaxies: intergalactic medium - QSOs: Q 0302-003 
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1. INTRODUCTION 

Once-ionized helium is the most abundant absorbing ion in the intergalactic medium 
(IGM). It outnumbers H I by a factor equal to the ratio of the ion number densities 
n-Reii/nm ^ 100, and therefore serves as an ideal tracer of the IGM in the early universe 
{z > 2). The presence of He II is signaled by absorption of He II Lya A304 A redshifted 
to the far-ultraviolet (far-UV) where it can be observed by space observatories such as 
HST. If He II ions reside in clumps or clouds positioned along the line of sight to a quasar, 
they produce discrete absorption lines similar to the H I Lya forest. But if instead they 
are diffused throughout the IGM, they will smoothly depress the flux level of a quasar 
shortward of its He II Ly a emission line. This is known as the Gunn-Peterson effect 
( punn fc Peterson 1965| ; [Scheuer 1965| ) . Although the Gunn-Peterson opacity was originally 



formulated for H I, there is an analogous formula for He II: 

, , _ C ?^HeIl(^) 0"HeII n 

where crneii = (vre^/meC^)/A is the resonant scattering Lya cross-section, and / and 
A are the oscillator strength and wavelength of the He II Lya line, respectively. The 
cross-section for H I is a factor of 4 larger than for He II. Hence, the ratio of opacities, 
R = Tueiiiz) / Tiii{z) = ri/4:, where r] = A^hcIi/A^hi, and A'hi and A'hcII are the column 
densities of H I and He II, respectively. In evaluating the Gunn-Peterson optical depth and 
throughout this paper, we assume Ho = 65 km s~^ Mpc^^, ^lo = 0.2 and A = 0. 

He II Gunn-Peterson absorption has been explored in a variety of numerical simulations 
(e.g. [Zheng fc Davidsen 1995| ; [Miralda-Escude et ai. I^ ; [Urott et al. TWTj [Zhang et al 



[19981; [Miralda-Escude, Haehnelt fc Rees 19991 , hereafter MHR). It has also been observed 



in four lines of sight. He II Gunn-Peterson absorption was first detected in the spectrum of 
the z = 3.286 quasar, Q 0302-003, obtained with the Faint Object Camera (Jakobsen et al. 
1994). However, because of the low spectral resolution and signal-to- noise ratio, the FOG 
data could not distinguish whether the absorption is produced by a uniformly distributed 
medium or by the He II counterpart of a forest of Lya lines (cf. Songaila et al. 1995). 
Subsequent spectra covering the wavelength range, 1140-1530 A taken by the Goddard 
High Resolution Spectrograph (Hogan et al. 1997) confirmed the strongly depressed flux 
level shortward of He II Lya in the QSO rest frame. It also showed a ledge in the flux level 
between 1280 - 1295 A, attributed to the proximity effect ( [Bajthk et al. 1988| ; Zheng & 



Davidsen 1995| ; [Giroux et al. 1995|) , an effect of increased ionization of gas in the vicinity 



of the QSO. In the region from 1240 to 1280 A, the GHRS spectrum was consistent with a 
null flux - however one with unacceptably large uncertainties (cf. Heap 1997). 



Measurements of the He II Gunn-Peterson effect provide a test of cosmological models. 
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Cosmological simulations can make concrete predictions of the epoch of IGM re-ionization 
from singly to doubly ionized helium (e.g. MHR), the timescale for the ionization change, 
and the characteristics of transition regions which might produce patchy He II absorption 
( [Reimers et al. 1997| ). Combining observations and simulations would constrain the density 
and ionization state of the gas, and the diffuse far-UV background radiation at z ~ 3. 

In this paper, we present new observations of Q 0302-003 made with the Space 
Telescope Imaging Spectrograph (STIS). STIS brings several improvements over GHRS 
for faint, point-source spectroscopy: first, a lower instrumental background resulting in 
higher S/N ratios; second, an imaging format that allows direct measurement of the (sky + 
instrumental) background; and finally, a wider spectral range covered in a single exposure, 
presenting an opportunity to study the properties of the IGM over a broader range in 
redshift. The STIS observations presented here permit the strictest constraint to date for 
He II Gunn-Peterson absorption, over a wide spectral interval, and thus offer a significant 
advance in our understanding of the physical processes at z ~ 3. In §0 and §|^, we describe 
the new observations. We then analyze in detail three important aspects of Gunn-Peterson 
absorption: the opacity of the absorption trough (§^, the "proximity effect" (§|^), and 
the appearance of opacity gaps in an otherwise solid wall of absorption (§|D. In §|^, we 
summarize the implications of our findings. 



2. OBSERVATIONS AND REDUCTIONS 

Our analysis of He II Gunn-Peterson absorption along the line of sight to Q 0302-003 
relies on two sets of observations: a far-UV spectrum of the QSO taken with STIS at 
low resolution, and high-resolution Keck spectra obtained by Songaila and by Ranch and 
Sargent, which were kindly made available to us. In this section, we describe the two 
datasets and our techniques of data reduction. 



2.1. STIS UV Spectra 

We obtained STIS observations of Q 0302-003 in December 1997 during the course 
of two "visits," each five orbits long (Program 7575). See Woodgate et al. (1998) and 
Kimble et al. (1998) for details of the instrument and its performance. The observation 
used the G140L grating which produces a spectrum covering the wavelength interval, 
1145-1720 A, at a nominal two-pixel resolution of 1.2 A. The sensitivity decreases strongly 
with wavelength below its peak at 1300 A: it is only 10% of its peak value at 1150 A, less 
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than 4% at 1145 A, and zero at at 1140 A. The total exposure time (10 orbits) was 23,280 s. 

Because of the faintness of the QSO, we selected a rather wide entrance slit (0'.'2) in 
order to maximize throughput. As a consequence, the resolution is degraded to about 
3 pixels, or 1.8 A. A wide slit also has the disadvantage of transmitting rather strong 
geo-coronal emission of H I Lya and a triplet of O I centered at A1302 (cf. Figure 1). 
This can be a serious problem because Lya lies in the Gunn-Peterson trough, and O I 
A1302 lies at its edge and might therefore interfere with the assessment of the "proximity 
effect". To counteract this problem, we obtained the spectra in the "TIME- TAG" mode 
with the intention of using only those data recorded during the earth-shadow time. We 
found, however, that contamination by geo-coronal emission can be fully accounted for in 
the reduction. We therefore assembled all of the TIME- TAG data from each orbit into a 
conventional spectrogram. 

The observations were reduced at the Goddard Space Flight Center with the STIS 
Investigation Definition Team (IDT) version of CALSTIS (Lindler 1998), which allowed 
us to make a customized treatment of the background. Such flexibility is needed in 
order to ensure accurate fluxes and opacities in the Gunn-Peterson trough. As shown 
in Figure |T], the background in a STIS G140L spectrogram is highly non-uniform. The 
strong geo-coronal emission at Lya and O I A 1302 affects a region wider than the projected 
entrance slit because of grating scatter. An additional source of background is the bright, 
diffuse background at the short- wavelength region of the image (the MAMAl "hotspot" ; cf. 
Landsman 1998). The QSO spectrum intersects this hotspot roughly through the middle. 
Because of the non-uniformity of the background, we sampled the background in two zones 
on either side of the QSO spectrum as closely as possible to the QSO (30-pixel offset) using 
an extraction slit height 30 pixels high (red lines in Figure |I|). In order to avoid smearing 
out the geo-coronal emission lines, we smoothed the background (15-point running mean, 
executed twice) only in the regions away from the emission lines. In contrast, the STScI 
data processing system (RSDP) samples the background of STIS G140L spectra within 
bands located ±300 pixels away from the QSO spectrum (inside each pair of dashed red 
lines in Figure |l]) where the background is only half that of our measured background. The 
resulting spectrum shows a spurious residual flux in the Gunn-Peterson absorption trough. 

Figure ^ shows our reduced spectrum overplotted with an error spectrum that assumes 
a/ZV statistics. Despite the low S/N ratio at short wavelengths, the residual flux below 
1175 A is real, since we see repeated appearances of signals at about the 5a level. There is 
also a feature just shortward of geocoronal Lya detected at the 5a level. The wavelength 
uncertainty is much harder to estimate. We made an on-board calibration between 
successive exposures, which ensures that the dispersion is precise. However, there is still the 
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Fig. 1. — STIS observations of Q 0302-003. The top panel shows a major portion of the 
co-added spectrogram (hnes 310-970 on the detector format). The two bright emission hnes 
are geo-coronal Lyman a (left) and 01 A1302 (right). The spectrum of the QSO has a strong 
continuum flux at A 1280 A, while the Gunn-Peterson absorption trough sets in at shorter 
wavelengths. Within this trough, there is a major flux flnger at A1231, and several minor 
regions of flux are visible shortward of Lyman a. The dashed red lines delimit the regions 
used to determine the background in the STScI reduction (not used for this paper), while 
the sohd red hnes apply to the region we adopted. The bottom panel shows the extracted 
spectrum, including the gross spectrum (QSO + background) shown by the black line and the 
average of the background above and below the quasar (squiggly red line). In the reduction, 
the background was smoothed (smooth orange line). Since the STIS MAMAl detector is a 
photon-counting imaging detector, the errors in the total counts are simply the square root 
of the counts. 
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possibility of a shift in the spectral direction between the QSO spectrum and the spectrum 
of the cahbration source. Such a shift might be induced by imperfect centering of the QSO 
within the 0'.'2 sht (8 pixels wide) or by thermally-induced drifts between the QSO and 
calibration exposures. Indeed, we find a 3-pixel shift in the zeropoint of the wavelength 
scale between the first exposure in a visit and the other four. Because of these uncertainties, 
we allow for a ±1 pixel (0'.'025) offset in the observed wavelength scale. Consistency with 
the wavelength scale of the H I Lyman forest suggests that a 1.06- A shift (1.8 pixels) is in 
order. 

The spectral resolution is also difficult to estimate because there are no sharp, well 
defined lines in the co-added spectrum. We therefore used the cross-dispersion profile of 
the co-added QSO spectrum as a proxy for the line spread function. The cross-dispersion 
profile grows broader toward shorter wavelengths. We estimate that in the Gunn-Peterson 
absorption trough, the FWHM is 3.1 pixels; thus the corresponding spectral resolution is 
1.80 A. 



2.2. Keck HIRES Spectra 

We used spectra of Q 0302-003 taken with the Keck HIRES spectrograph (Vogt et 
al. 1994) provided by A. Songaila and by M. Ranch and W. Sargent to compare H I Lya 
absorption with the He II Lya absorption in the STIS spectrum. Table |l| summarizes the 
data. Successive columns give the observed wavelength range of the echellogram in A, the 
equivalent wavelength range apphcable to the STIS spectrum (i.e. the observed wavelength 
divided by 4 for the He II /H I wavelength ratio), the observed resolution in km s~^, 
the maximum signal-to-noise ratio, and observer /source of the data. Hu et al. (1995) 
and Ranch et al. (1999) provide reduction details for the Songaila and Ranch datasets 
respectively. We processed the Ranch data and estimated the la uncertainties from the 
detector noise and photon counting statistics; for the Songaila data we calculated the la 
uncertainties from the apparent rms noise in the spectrum itself. In our analysis, we relied 
on the Ranch spectrum at A ^ 4600 A where it has higher S/N ratio; we used the Songaila 
spectrum at shorter wavelengths for some metal-line identifications and to improve our H I 
Lya profile fits using higher members of H I Lyman series. We detected a 2-3% residual 
fiux in the bottoms of saturated absorption lines in the Songaila data, but found that it 
makes no significant difference to our analysis. 
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Fig. 2. — The STIS spectrum of Q 0302-003 after applying a background correction. Top: 
the absolute flux distribution in erg s^^ cm~^ (solid line) showing He II absorption 
shortward of 1300 A, overlaid with the 1 a uncertainty (dashed line). Bottom: a magnified 
portion of the spectrum in the Gunn-Peterson absorption trough. In this panel, the flux was 
normalized to unity at — 2.1 x 10~^^ erg s~^ cm~^ A~^, and binned by 3 pixels in order 
to show the average flux in a resolution element. 
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Table 1: Keck HIRES Spectra of the H I Lyman Forest Toward Q 0302-003 



^vac 


(A) 


A/4 


R (km/s) 


(5/iV)^ax 


Observer 


3650- 


-5983 


912-1496 


8.3 


70 


Songaila 


4220- 


-6656 


1050-1664 


6.9 


130 


Rauch & Sargent 



Absorption lines of H I, C IV and other ions detected in the Keck spectra were 
fitted with Voigt profiles with the program, VPFIT (Webb 1987). We base our redshifts 
on wavelengths corrected to the vacuum heliocentric frame. Higher members of the H I 
Lyman series were used wherever possible to constrain the profile fits. We estimate the 
completeness limit at logiVni ~ 13.0 (cm~^) based on the point at which the H I column 
density distribution turns over. Our results for redshift Doppler parameter 6hi, and H I 
column density, A^hi, largely agree with those published by Hu et al. (1995). The few 
differences that we find can be attributed to {€) the higher resolution of the Rauch dataset, 
[ii) newly identified metal lines, or {Hi) an accounting for the ink spot in the middle of 
the HIRES CCD. The absorption-line parameters for the H I systems will be used in §^ to 
model the He II absorption. A detailed line list will be presented elsewhere ( |Rauch et al. 
19991) . 

3. OBSERVATIONAL RESULTS 

In this section, we present the STIS spectrum of Q 0302-003, which covers the 
wavelength interval, 267-402 A in the rest frame of the QSO. We make this presentation 
in the form of two comparisons. In §3.1, we compare the He II Lya absorption trough 
(Arest < 304 A) with the corresponding H I Lya forest. In §3.2, we compare the He II Lya 
absorption in the spectrum of Q 0302-003 with other lines of sight. 

3.1. Comparison of the He II and H I Lya Spectra 

In order to measure the strength of He II Lya absorption, it is necessary first to 
define the continuum fiux distribution of the QSO in this spectral region. Of course, the 
continuum is not a directly observable quantity, since it has been blocked by He II Lya 
absorption. We must therefore extrapolate the fiux longward of the Gunn-Peterson trough. 
But even this continuum is contaminated by absorption by intervening systems. Figure § 
shows a "widened spectrum" of the QSO in order to highlight the absorption features 
longward of the Gunn-Peterson trough. (It was constructed by vertically shifting successive 
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raw spectrograms by 2 pixels and superposing them.) There is a break in the spectrum 
at ~1620 A suggestive of a Lyman hmit system at z = 0.773. We searched for its Mg II 
A2800 absorption but could not conclusively identify these metal lines, because they are 
obscured by strong H I lines. There is also a known metal-line system at z = 1.889 (Hu et 
al. 1995). A strong absorption line {W\ = 1.2 A) at 1689 A is probably He I Lya A584 
from this system. Because of the absorption from these intervening systems, the QSO 
flux distribution is not a smooth power law and hence, not easily extrapolated to shorter 
wavelengths. 

The lower panel of Figure |^ shows the observed flux distribution with two continuum 
estimates superposed. One attempts to account for the absorption of intervening 
systems via the program, CLOUDSPEC (Hubeny 1998), a combination of CLOUDY 
( l^erland et al. 1998| ) and SYNSPEC (Hubeny et al. 1994). The other (and the one 



we adopt) is a flat continuum in wavelength (corresponding to a = 2, ~ z/ ") set at 
Fc=2.1 X 10^^^ erg s~i cm'^ A~\ 

Figure ^ compares the STIS spectrum of Q 0302-003 with a simple model of the He II 
Lya forest derived from the high- resolution (6 km s~^) Keck/HIRES spectrum of the H I 
Lya forest. The model calculation involved four steps: (i) estimate the optical depth in 
the H I Lya lines from the flux of each data-point normalized to the apparent continuum, 
thi(A) = — ln(/™); (ii) calculate the corresponding optical depth in the He II lines rHeii(A) 
at each data-point assuming a ratio of He II to H I optical depths R = THeii(A)/rHi(A); 
(iii) reconstitute the He II Lya forest, /(A/4) = exp(— rHeii(A)); and (iv) degrade the 
spectrum to the STIS resolution (1.8 A). The values of i? = 25 and 100 used in the figure 
correspond to a ratio of He II to H I column densities, i] = Nueii/Nm = 100 and 400 
respectively, and assume pure turbulent line broadening. The comparison makes clear what 
strong absorptions are part of the Lya forest and what absorption is associated with the 
diffuse, or underdense, medium. For example, the high opacity of He II Lya in the region, 
1260-1280 A, known as the Dobrzycki-Bechtold void (pobrzycki fc Bechtold 1991|; hereafter 



the D-B void), must originate in the diffuse medium. On the other hand, the H I opacity 
gap at 1231 A is also a gap in He II. At the very shortest wavelengths, there is a fair 
correlation of the observed spectrum with the simulated He II Lya forest, suggesting that 
the opacity of the diffuse medium may be lower in this region. Close to the QSO redshift, 
there is also significant flux that has been interpreted as evidence of the proximity effect 
(Hogan et al. 1997). 
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Fig. 3.— Far-UV flux distribution of Q 0302-003. Top: a 2-D display of the 10 spectra 
stacked verticaUy to show the reaUty of spectral features. Bottom: the observed flux 
distribution in units of erg s~^ cm~^ A~^. Two contimiTim estimates are shown as gray 
lines: one allowing for a Lyman limit system at z = 0.773, and the other set to a constant 
value. The latter was used to normalize the flux in the Gunn-Peterson trough. 
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Fig. 4. — Comparison of the observed (black line) and simulated (filled gray) spectra 
of Q 0302-003 in the Gunn-Peterson trough. The simulations assumed pure turbulent 
broadening (^ — 1) and two different values of the He II/H I optical depth ratio, R = 
25 and 100. The region contaminated by geo-coronal Lya was omitted in order to avoid 
confusion. 
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3.2. Comparison With Other Lines of Sight 

To put the line of sight to Q 0302-003 in context, Figure |^ compares the STIS spectrum 
of Q 0302-003 with those of PKS 1935-692, a z = 3.18 quasar observed by the GHRS 
(Tytler et al. 1995, Jakobsen 1996) and by STIS (Anderson et al. 1999), and HE 2347-4342, 
a. z = 2.885 QSO observed by GHRS (Reimers et al. 1997). (The raw observational data of 
PKS 1935-692 and HE 2347-4342 were re-reduced at Goddard.) The abscissa in this figure 
is now the redshift of He II Lya, z = (A/303.78) — 1, instead of observed wavelength. In 
each panel, the downward pointing arrow shows the wavelength of He II Lya at the redshift 
of the QSO. Below, we discuss the three major features of the He II Lya absorption trough: 
the proximity effect, the absorption trough itself, and gaps in the absorption. 

The Proximity Effect. This term refers to an observed decrease in the number of 
H I Lya - absorbing clouds in the neighborhood of a QSO, at least in part caused 
by the H-ionizing radiation field of the QSO. A He II proximity effect involving both 
He Il-absorbing clouds and diffuse components of the IGM is also expected. It is best 
defined by the spectrum of Q 0302-003, although it is also evident in PKS 1935-692; 
however, note that the edge of the proximity zone is blocked by a pronounced gap in the 
absorption. In §|], we will use the observed fiux distribution in the proximity zone to explore 
the shape of the radiation field of the QSO and the metagalactic background at z ^ 3. 

The He II Gunn-Peterson Absorption Trough. Strong, continuous absorption is evident 
in all three quasar spectra outside of the proximity zone. Of the three, only the spectrum 
of Q 0302-003 has the broad baseline in redshift [z = 2.78 — 3.28) needed for studying 
the evolution of the IGM. The spectrum of PKS 1935-692 would appear at first also to 
be useful for this purpose, but it is contaminated by a z = 0.30, damped Lya system 
(logA^Hi = 21.2) whose Lyman continuum absorption effectively blocks the light from the 
quasar corresponding to He II Lya at z < 3.0 (the shaded region in Figure 5). In §4, we will 
analyze the absorption trough in the spectrum of Q 0302-003 in detail. 

Gaps in the Absorption Trough. At wavelengths shortward of the shelf of transmitted 
fiux that we attribute to the proximity effect, there are significant gaps in an otherwise 
solid wall of absorption by He II. As can be seen in Figure ^ there are windows of 
transmitted fiux in the spectrum of Q 0302-003 at z = 3.05, PKS 1935-692 at 3.10, and 
HE 2347-4342 at z = 2.87 and 2.82. Since these windows of transmitted fiux are broader 
than our instrumental resolution and are not accompanied by similar, much narrower 
ones elsewhere, it is doubtful that they correspond to gaps appearing by chance among 
randomly distributed absorbing clouds. In we shall argue that these coherent stretches 
of transmission are probably the result of discrete sources photo-ionizing the surrounding 
material. 
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Fig. 5. — Comparison of the He II Gunn-Peterson absorption trough plotted as a function 
of redshift for the He II Une at 303.78 A. Q 0302-003 and PKS 1935-692 were observed 
with STIS, and HE 2347-4342 by GHRS. The downward arrows denote the redshift of the 
QSO. As in Figure 4, the spectral region near geo-coronal Lya is grayed out. In the case 
of PKS 1935-692, which has an intervening damped Lya system a,t z — 0.30, the predicted 
absorption spectrum from the DLA is also shown in gray as a warning that the opacity at 
z < 3 should not be interpreted as Gunn-Peterson absorption. 
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Note how the morphology of low-opacity regions changes with redshift. At redshifts 
z > 3 (and outside of the proximity zone), the opacity gaps are isolated, discrete, and major. 
In contrast, regions of transmission at 2; < 2.9 occur more frequently but with smaller 
amplitude. This changing morphology will be discussed in the following three sections. 

4. THE He II GUNN-PETERSON ABSORPTION TROUGH 

In this section, we make quantitative assessments of the He II Lya optical depth of 
the IGM along the line of sight to Q 0302-003. These include both direct determinations 
of the average transmission within the He II Gunn-Peterson absorption trough (§4.1), and 
a more sensitive but less direct method that makes use of the low level of fluctuations in 
transmission in the trough (§4.2). Since the measurements cover a wide redshift range 
{z = 2.78 — 3.28), we were able to trace the redshift evolution of the He II Lya opacity. 
We find that it rises rapidly with redshift, a result that is compatible with a shift in the 
C IV/Si IV ratio between z = 2.9 and z = 3.0 reported by Songaila (1998) (§4.3). On the 
assumption that a lower He II Lya opacity is due to a higher fraction of doubly-ionized 
helium, we follow the re-ionization history of helium over the redshift interval, z = 2.8 — 3.2 
(§4.4). 

4.1. Measured Opacities 

We estimated the He II Lya opacity for different regions as shown in Figure ^ The 
upper panel shows the normalized flux (i.e. the relative transmission of the IGM) as a 
function of redshift, while the bottom panel shows the spectrum of apparent optical depths. 
The region labeled "D-B" is the Dobrzycki-Bechtold void (1991). Regions A and B are 
used together to characterize the He II opacity at 2; ~ 3 for comparison with theoretical 
predictions. The region labeled "S" corresponds to redshifts that Songaila (1998) identified 
as the high-ionization regime (cf. §4.3). We measured the average residual intensities in 
each region (/) and converted them to representative optical depths (r = — In J) with results 
listed in Table ^ Successive columns give the region ID, the number of data-points in the 
region, the average transmission and its uncertainty, the rms error of a single data-point, 
and the optical depth corresponding to the average of the transmitted intensities. 
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Fig. 6. — Transmissions (upper panel) and opacities (lower panel) in the Gunn-Peterson 
absorption trough. The selection of the different spectral regions is explained in the text. 
The measured values of the apparent transmissions and optical depths are given in Table H 
and depicted by gray bars. 
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Table 2: Measured Opacities in the Gunn-Peterson Trough 



Region 


z 


N 


I 




^region 


D-B 


3.18 


14 


0.0114 ± 0.0044 


0.0164 


A 47+0.48 


A 


3.15 


73 


0.0086 ± 0.0029 


0.0246 


^■''-'-0.29 


B 


2.93 


52 


0.0151 ± 0.0079 


0.0567 


4.i9i°:l 


S 


2.82 


52 


0.1532 ± 0.0179 


0.1292 


1.881°:}? 


A+B 


3.06 


125 


0.0113 ± 0.0037 


0.0410 


4.48l°:i 



In agreement with previous GHRS observations, the optical depth of the D-B void 
Tdb = 4.471q'33 is quite high, implying that the absorption appears to be due to diffuse or 
underdense material. Alternatively, there may be a nearby starburst galaxy or other soft 
ionizing source that could ionize H I but not He II, thereby producing a low H I opacity but 
not affecting the He II opacity. As the GHRS spectrum only sampled the Gunn-Peterson 
trough above z — 3.1, the measured opacities in the other regions (A, B, S) are new, and 
will be discussed later in this section. 

The rms errors quoted in the table give a measure of the random errors. There 
remains the question of systematic errors, which are the more serious of the two. Since the 
STIS detector is a pulse-counting detector with a hnear response at low light levels, we can 
rule out the possibility that the low residual flux in the Gunn-Peterson trough is due to a 
systematic error in the response function. However, there could be a systematic error in the 
residual flux if the cross-dispersion ( "vertical" ) gradation in the background were not linear, 
so that the average of the upper and lower background regions would not be representative 
of the background under the QSO spectrum. We can put an upper limit on this systematic 
error by re-reducing the spectrograms with the upper background or lower background only. 
These alternate reductions produce residual fluxes in the Gunn-Peterson absorption trough 
that are well within the uncertainties quoted in the table. Other reduction procedures 
(15-point background smoothing applied twice, registration of the observations before 
co-addition) yield differences on the same order as adopting the upper or lower background 
only. We therefore interpret the errors listed in the table as representing total errors. 

4.2. Opacity Fluctuations 

Now that we have explicitly measured the opacities, we can obtain a more sensitive, 
but less direct result by looking for gaps in the opacity that appear by chance and 
comparing them with predictions. We find that a high opacity is required to produce the 
broad stretches of nearly zero transmission as seen, for example, in spectral region A. This 
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region has 24 independent samples of width iS.z = 0.0058, where independent samples are 
assumed to be separated by the instrumental resolution (3 pixels). All 24 samples show a 
relative transmission, I < 0.030 (cf. Figure 2), implying that the optical depths exceed 
— ln(0.03) = 3.50 everywhere in this region. We designate this opacity as rumit- Even in 
the absence of extra ionizing sources, we should expect to see some variations caused by 
the random gaps in the overlapping of clouds. The fact that such variations do not cause 
observable transmission spikes implies that the average optical depth is greater than a 
certain amount - a value that we shall now evaluate. 

Fardal, Giroux & Shull (1998) constructed a model for the clouds in the IGM based on 
the observed properties of the Lya forest and then, on the basis of Monte Carlo simulations, 
found that fluctuations in the apparent opacity of He II, defined as r = — In /, over an 
interval Az, follow the approximate relationship (their Eq. 17): 



T / model V 4 / 



where f is the average of r over many intervals of Az, and At represents a deviation of r 
from this average. The quantity ^ is the characteristic ratio of the He II line width to that 
of H I (ranges from 0.5 to 1.0, depending on whether the line broadening is dominated by 
thermal or turbulent processes), and ?7ioo = A^Heii/(100 Nm). 

Now consider the (very small) probability p that a single sample of width Az shows a 
spike I > 0.030, i.e. rumit < 3.50. This quantity p then represents the one-tail area below 
Tiimit for a normal distribution with a standard deviation At centered on f . The chances of 
not seeing a spike over all 24 samples is e~^^ ^. We now determine the limiting condition 
that would create a 10% expectation of seeing no spikes in any of the samples, since this 
leads to a confidence level of 90% that we are not overestimating f and/or underestimating 
At. This happens when p = — ln(0.1)/24 = 0.096 which is the one-tail area for a = 1.30 
standard deviations below the center of the normal distribution. This requirement sets a 
lower limit for the average opacity. 



T > Tlimit 



1 



T / model 



-, -1 



(3) 



For Tiimit = 3.50, cr = 1.30, ?7ioo =3.5 (cf. § |5.4.3| ), and ^ = 0.5 or 1.0, we obtain r > 6.6 and 
6.2, respectively. Both values are somewhat higher than the measured value, ta = 4.75 ±0.4 
in Region A. 
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4.3. A Break in He II Lya Opacity at z = 2.9 - 3.0? 

Songaila (1998) recently reported on a study of Keck HIRES spectra of 13 quasars, in 
which she found that the C IV/Si IV ratio aX z = 2.9 was 3.4 times higher than aX z = 3.0. 
Since the ionization potentials of Si IV (45.1 eV) and C IV (64.5 eV) straddle the ionization 
potential of He II (54.4 eV), C IV/Si IV serves as a proxy for the ionization state of helium, 
and she attributed the change in this ratio to a hardening of the metagalactic spectrum 
below z = 3.0. However, Songaila's findings were disputed by Boksenberg et al. (1998) 
who found no sudden change in ionization, although they did note a gradual increase in 
C IV/Si IV toward lower redshifts between z ^ 3.8 and z ^ 2.2. 

The STIS spectrum of Q 0302-003 provides an opportunity to trace the ionization 
evolution of helium over a broad redshift range, z = 2.8 — 3.2. (We assume that a lower 
He II opacity means a higher fraction of doubly ionized helium.) Our calculations with 
CLOUDY ( [Ferland et al. 1998|) indicate that the He Ill/He II ratio should change by at 
least as much as C IV/Si IV for a Madau et al. (1999) background radiation field at z = 3. 
We therefore looked to see whether the He II opacities show a distinct change between 
z = 2.9 and z = 3.0. We find that Region S of Figure ^ does indeed show a residual flux in 
the Gunn-Peterson trough at z < 2.87 implying an average transmission that corresponds 
to Ts = 1.88 in this region. 

The redshift interval, z = 2.77 — 2.87, however, is by no means a void in the 
conventional sense. Many H I Lya forest lines occur there, in some cases with clear 
corresponding increases in He II opacity. To compare the z < 2.9 and z > 3.0 opacities, we 
selected those data points in Region S unassociated with features in the Lyman forest, i.e., 
positions identified as having a normalized H I transmitted flux greater than 0.95 even when 
degraded to STIS resolution. We find 10 non-contiguous data-points in Region S fulfilling 
this /hi > 0.95 requirement. The He II opacity corresponding to the average transmission 
for these 10 data-points is r = 1.41 ± 0.06, which is a factor of 3 times lower than in the 
D-B void. Similarly, the opacity for the whole Songaila region {z = 2.82) which includes 
numerous Lya lines, Ts = 1.88, is a factor of 2.5 times lower than for Region A at z = 3.15. 
We conclude that the He II Gunn-Peterson data are in accord with Songaila's assertion of 
an abrupt shift in ionization between z = 2.9 and z = 3.0. 

4.4. Comparison with Cosmological Models 

Observations of He II Lya absorption give a direct connection with cosmological 
simulations if differences in resolution are taken into account. For example, simulations by 
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Zhang et al. (1998) cover a cube only 9.6 comoving Mpc on a side with a cell size of 75 
kpc^. In contrast, our STIS spectrum of the G-P trough covers a much larger distance (400 
comoving Mpc) but at a coarser resolution (~ 2 comoving Mpc), so it smooths over all of 
the fine-scale fluctuations predicted by theory. Nevertheless, we can apply the same kind of 
smoothing on the theoretical opacities by taking averages. Zhang et al. (1998, their figure 
20) predict that the average He II opacity should increase from f = 0.50 at z = 2.4, to 
f = 0.83 at ;z = 3, and f = 2.43 at z = 4. 

The measurements listed in Table 2 imply that the total He II opacity rises more 
rapidly with redshift than the models (Zhang et al. 1998, Fardal et al. 1998). As shown in 
Figure |^ which combines data for HS 1700+6416 ( [Davidsen et al. 1996| ) with the measured 



opacities derived from Regions S, B, and A in the spectrum of Q 0302-003, the He II optical 
depth in Lya increases sharply from r = 1.0 at z = 2.4 to r ~ 4.7 at 2; = 3.15. The more 
sensitive estimates of the He II opacity based on the lack of opacity fluctuations suggest 
that it may rise with redshift even faster than shown. The top panel also plots the mean 
He II Lya opacity predicted by Zhang et al.'s (1998) cosmological model. It shows that the 
model generally underestimates the He II Lya opacity and does not reproduce the observed 
steep rise in opacity with increasing redshift. 

The top panel of Fig. |^ also shows the opacity expected if all helium were in the 
form of He"^ and were uniformly distributed. Two curves are shown: one for our assumed 
cosmological parameters (solid line), the other for Madau et al.'s (1999), i.e. Ho=50 
km s~^Mpc~^, Qq = 1.0 and A = 0. In either case, the measured opacity is more than 
600 times smaller than expected. The fact that we detect a small residual intensity in 
Regions A and B indicates that most, but not all, of the helium is doubly ionized at z = 3. 
Accordingly, the bottom panel shows the implied fractions of He III and He II. 

In summary, away from distinct gaps in the He II absorption, the Lja opacity shows a 
sharp rise between z = 2.8 and z = 3.2. At 2; ~ 3.15, the optical depth is estimated at 4.7 or 
more. It therefore appears that we are witnessing the transition from a He II Lya forest to a 
He II Gunn-Peterson trough. The sharp rise in opacity at 2 ~ 3 (interpreted as an increased 
fraction of He II ) is consistent with Songaila's claim of a lower ionization level of the IGM 
at z > 3 . However, we caution that these results are based on a single line of sight. Other 
lines of sight may give different results since He II is such a minor fraction of helium (0.16%). 
In fact, some of the large scatter in Songaila's measured C IV/Si IV ratios may reflect such 
differences. We plan to pursue the question of a "universal epoch of reionization" by studies 
of the He II Gunn-Peterson effect along the sightlines to PKS 1935-692 and HE 2347-4342. 
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Fig. 7. — Evolution of the He II Lya opacity (top) and helium ionization fractions (bottom). 
Top: The filled circles represent data from our STIS observations of Q 0302-003, while the 
filled triangle represents the HUT results for HS 1700+6416 (Davidsen et al. 1996). The 
bold curve shows the opacity predicted by Zhang et al. (1998). The other curves show the 
expected He II Lya opacity if all helium were singly ionized and uniformly distributed. The 
assumed cosmological parameters are: Ho=Q5 km s^^Mpc^^ and Oq = 0.2, A = (solid 
line) , or Ho=50 km s~^Mpc~^ and Qo = 1.0, A = (dash-dot). Bottom: Implied helium 
ionization fractions in the redshift range observed by STIS, based on the assumption that 
all hehum is once or doubly ionized. The fraction of Hell is only about 0.1 % at 2; ~ 3. 
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5. MODELS OF GUNN-PETERSON ABSORPTION AND THE 

PROXIMITY EFFECT 

The He II opacities measured in §^ imply that the vast majority of He is doubly ionized 
at 2; ~ 3, presumably by the UV background. In this section, we derive the properties 
of the UV background at 2; ~ 3 from the H I and He II Lya spectra of Q 0302-003. To 
pursue the matter, we constructed spectral models that take into account just one ionizing 
source, the QSO itself, in addition to the UV background. These models allow us (1) to 
identify regions within the proximity zone that are particularly sensitive to the shape of 
the UV background, and can thus be used to constrain it, and (2) to discriminate among 
the He II opacity gaps those likely caused by underdense regions in the IGM from those 
due to discrete ionizing sources. In the following sections, we describe the model input, 
calculations, and results. 

5.1. Model Input 

The three main inputs to the models are the observed flux distribution in the He II 
Gunn-Peterson absorption trough, the QSO flux distribution, and the physical parameters 
of the H I Lya forest. Below, we describe our methods for estimating the QSO continuum 
flux distribution and the physical characteristics of the H I Lya lines. 

5.1.1. The spectral energy distribution of Q 0302-003 

We estimated the ionizing flux distribution of the QSO from measurements of the 
apparent continuum combined with model continuum flux distributions. Figure |^ plots 
the theoretical luminosity distribution of Q 0302-003 along with observationally derived 
continuum points transformed to the quasar rest frame [z = 3.29). The observed points 
are taken from our STIS data and from the ground-based spectrum obtained by Sargent, 
Steidel, & Boksenberg (1989). The STIS observations cover the restframe interval, 267 - 
402 A, or 7.5 x 10^^ - 1.12 x 10^^ Hz. The SSB fluxes were muhiplied by a factor of 1.25 to 
account for the depression of the apparent continuum by the multitude of lines that make 
up the Lya forest. This factor was estimated from the normalized Songaila et al. spectrum 
(Table 1) degraded to the resolution of the spectrum obtained by Sargent et al. We found 
that a power law of the form /Q^o _ j^QSO ^ ^^^-^-a^^ with a^f^ = 1.9, gives a good 
representation of the data from the optical to the far-UV (Arest = 1000 — 304 A). In our 
models of the QSO flux distribution, we extrapolate this power law into the unobserved 
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He II Lyman continuum as well. 

The two theoretical models shown in the figure are for a bare accretion disk around 
a super-massive, Kerr black hole, one with a black- hole mass M = 2 x 10^° Mq, and 
accretion rate M = 24 MqJt~^ (upper curve); and the other with M = 1.6 x 10^*^ Mq and 
mass accretion rate M = 16 Moyr~^ (lower curve). Both models assume the maximum 
stable rotation of the black hole, with a specific angular momentum (a/M) = 0.998. The 
integrated spectrum of the disk taking into account all general relativistic effects was 
computed by program KERRTRANS ( |Agol 199^ ); the best fit was obtained for disk seen 



almost face-on. A detailed description of the modeling procedure is given by Hubeny & 
Hubeny (1998); the presented models are taken from the grid of Hubeny, Blaes & Agol 
(1998). As these spectra represent a bare accretion disk, they do not take into account the 
effects of a Comptonizing corona. 

The theoretical flux distribution has three distinct spectral regions, each with its 
characteristic spectral slope, a, deflned by / oc z/~". For frequencies z/ < 3 x 10^^ Hz 
{Xrest > 1000 A), a ^ 0.55. For z/ = 3 x 10^^ - 1 x 10^^ Hz {Xrest = 1000 - 300 A), a ^ 2. 
At z/ > 10^^ Hz, the flux falls off steeply, and a ~ 11. The models reproduce the observed 
flux in the flrst two regions. In the high-frequency region, where we have no observed data, 
the theoretical flux is likely a lower limit, because the models do not take into account 
the effects of a Comptonizing corona, which increase the flux dramatically. Therefore, we 
conclude that the power-law slope of the He II Lyman continuum in the QSO spectrum is 
2^<^°^11. 



5.1.2. Line lists 

In order to estimate the value oi rj = A^Heii/A^Hi as a function of z, we made use of 
a line list containing the observed Lya lines whose parameters (2;, A^hi, ^hi) were derived 
from VPFIT measurements of the Keck HIRES spectra (cf. § p.2|) . This line list is hereafter 
referred to as the 'observed line list'. Following Songaila et al. (1995), we expanded the 
sample in the range of 9 < log [Nm) < 13, by producing a line list with random values of 
Ahi; b and z, whose distribution parameters were obtained by interpolating in redshift the 
values determined by Kim et al. (1997). The resulting line list is called the 'combined line 
list'. Note that even ii r] = 1000, the log Ahi = 9 clouds do not contribute much to the He II 
opacity. 

The immediate effect of including this plethora of very weak lines is to depress the 
apparent continuum in the H I Lya forest, i.e. to produce a shallow H I Gunn-Peterson 
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Fig. 8. — Probable spectral energy distribution of Q 0302-003. The filled squares show 
the observed continuum fiux transformed to the quasar rest frame. The upward-pointing 
arrows give the frequencies of HI Lja (left) and the HI Lyman limit (right). The downward- 
pointing arrows show the same quantities for He II. The sohd lines show two theoretical 
models of a bare accretion disc around a super-massive, Kerr black hole with a black-hole 
mass M = 2 X 10^° Mq and the mass accretion rate M = 24 Mq yr~^ (upper curve); and 
M — 1.6 X 10^° Mq and mass accretion rate M = 16 M0yr~^ (lower curve). Both models 
assume an almost face-on disk, with inclination cosi = 0.99. 
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absorption trough. Figure ^ compares the H I Lya forest computed with the combined hne 
hst vs. the observed hne hst. The predicted Gunn- Peterson trough has a depth of about 
6%, which is shaUow enough to escape detection in real data. However, it is consistent with 
Fang et al.'s (1998) measurement (thi = 0.113 ± 0.02) of the H I Gunn-Peterson effect in 
the z = 3.787 QSO, PKS 1937-101. 

5.2. Model Calculations 

The model calculation is recursive. First, we consider the cloud closest to the QSO: 
we assume that it sees the 'pure' QSO spectrum and the UV background. We compute its 
H I and He II photo-ionization rates, whose ratio is proportional to the value of r] for this 
cloud (with the simplifying assumptions used), and allows the derivation of its He II column 
density. We then consider the second closest cloud to the QSO: it sees the UV background 
and the QSO whose spectrum shortward of the H I and He II Lyman limit is depressed by 
the continuum absorption due to the first cloud. Since we know the H I and He II column 
densities and redshift of the first cloud, we can compute the QSO spectrum seen by the 
second cloud, the resulting H I and He II photo-ionization rates; and from the observed H I 
column density of the second cloud, we compute its He II column density. The process is 
then repeated for each cloud of the considered line list. We then use the calculated value 
of A^Heii and the line-width ratio bueii/bm = 1 (pure turbulent broadening) to calculate the 
He II absorption spectrum (Voigt profiles) of all the lines in the line list. 

The main difference between the development made here and the theory of Stromgren 
spheres is the presence of a diffuse ionizing radiation field (the UV background) which 
leads to a much more gradual stratification in ionization compared to the sharp edge of 
a Stromgren sphere. Also, contrary to a Stromgren sphere, the gas is considered to be 
inhomogeneous, i.e. composed of numerous clouds embedded in a diffuse medium. Finally, 
there is a velocity gradient within the ionized region (due to the expansion of the Universe) 
so that clouds far from the main ionizing source, i.e. the QSO, recede faster from the QSO 
than do closer ones. 



5.2.1. Assumptions 

To compute the ratio of column densities of He II and H I, we assume that Lya clouds 
at 2; ~ 3 are in photo-ionization equilibrium with hydrogen mostly ionized and helium 
mostly doubly ionized. Indeed, since all the Lya clouds in front of Q 0302-003 have 
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Fig. 9. — Effect of very weak Lja lines on the observed H I Lya forest. Top: tlie normalized 
spectrum of the Lya forest computed with observationally derived parameters, A^hi, b and 
z. Bottom: the same as above except that the "combined" hne hst was used (see text). 
The effect of including very weak lines is to produce a shallow H I Gunn-Peterson trough by 
depressing the apparent continuum by about 6%. 
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logA/^Hi < 15.47, most H is certainly ionized, and Figure |^ shows that nearly all of the 
helium must be doubly ionized. We also assume that the ratio of column densities of He II 
and H I is given by 

-^Hell 'T-He ttHell ThI , 

iVHi nn ohi THeii 
where is the recombination coefficient of ion i (H I or He II), Fj is its ionization rate, and 
'^He, nn are the total number densities of hehum and hydrogen respectively. Note that for 
clouds with high A^Heii, self-shielding and emissivity of the clouds start to be substantial, 
and this simple equation is no longer valid. For rinc/nu = 0.082 and aHeii/«Hi = 5.418 (cf. 
Osterbrock 1974), the equation reduces to 

r, = OAA^. (5) 

J- Hell 

The photo-ionization rates are given by: 

F, = / a, du, (6) 

where / is the total (QSO + UV background) ionizing photon flux, cTj^ ~ Gi {ujui)''^ is the 
photo-ionization cross-section (cf. for example, Osterbrock 1974), and z/j is the frequency of 
ion i Lyman limit. At all redshifts, the photo-ionization rates due to the UV background, 
F^, are related to the intensities at the Lyman limits Jj by the relation: 

J ^ 47rai Jj 

' h 3 + ^' 

if the Lyman continuum radiation of the UV background spectrum can be represented by a 
power-law J(z/) = Ji{iJ / iJi)~l^i . Useful, numerical relations between these quantities are: 

ri, ^ 2.99 . 10-- ^ ^ 
ri„ - ^.47 X 10-.^ ^ ^ (8) 



5.2.2. The first cloud 

Starting with the closest cloud to the quasar and moving towards decreasing redshifts, 
we compute the expected value of rj as the contributions from the quasar decrease while the 
UV background remains the same. The photo-ionization rates of the first cloud aX z = zi 
are: 

h 3 + aP° I 1 + ^1 ) ' ^ ^ 
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where the QSO flux at a redshift z and frequency z/ > z/j is j'^'^'^i^z^v) = 
jQSO^^'j ^^-^ _|_ z//z/j]~"i . The quasar flux at the Lyman hmit seen by a 

cloud at redshift z is related to the QSO luminosity by: 



l + ^QSO 47rD£ 



where Di^{z,zqso) is the luminosity distance between the cloud and the QSO (cf. Kayser, 
Helbig & Schramm 1997). The photo- ionization rates of H I and He II are calculated with 
Eq. ^ and then Eq. || is used to compute rj. The value of the He II column density is then 
easily calculated from the observed H I column density. 



5.2.3. The kth cloud 

Since the QSO spectrum below the ion i Lyman limit seen by the kth cloud {k > 2) is 
depressed by the continuum opacity produced by the k — 1 clouds located between it and 
the QSO, the photo-ionization rates of the kth cloud at 2; = are: 



Ui/{1+Zi,) h V 



Oy exp 



E 



where t] = ai + z\)/vi)~'^. The values of N\^i are obtained from the observations 
and the values of iVfjeii have been estimated at the calculation of the l"^^ cloud. Expanding 
the exponential as a series (exp (— x) = l]m=o(~''')'"/'^0 obtain an analytical formula 
for the second term in Eq. O, which can be conveniently written as: 



Ji (^i n- + zqso 



h 



QSO 



E 



.ly 



ml (3m + + 3) 



k-i 

E 

.1=1 



N a,; 



1 + 

l + Zi 



(12) 



In practice, the number of terms to be evaluated depends on the number and values of the 
largest Nl lines; we find that a summation on m up to ~ 30 is sufficient for our purpose. 



5.2.4. Model spectrum 

Once the ionization rates are known, we can calculate rj (Eq. 5) and A^Heii- In the 
final step, we apply the calculated value of A^Heii and the line-width ratio bueii/bm = 1 
to obtain the He II absorption spectrum (Voigt profiles) of all the lines in the line list. 
After convolution by a line spread function of FWHM=3.1 pixels (cf. §2.1| ), we obtain a 
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model spectrum of the He II Lya absorption trough suitable for comparison with the STIS 
observations. By adjusting the input parameters in different simulations, we can measure 
diagnostic properties of the UV background and evaluate the need for diffuse gas in the 
IGM. 

To summarize, ri{z) can be evaluated based on only the observed values of A^hi? the 
QSO observed flux and spectral slope, and on some assumptions about the UV background. 
By matching the observed spectrum, we can set constraints on the H I photo-ionization rate 
due to the UV background, and on the softness parameter, 

S ^ (13) 

^ HcII 

We prefer to use S, which intrinsically takes into account the shape of the UV background 
spectrum, instead of another, often used measure of the shape of the UV background: 

= (14) 

o'Hell 

where Jhi and JhgII are the UV background intensities at the H I and He II Lyman limits, 
respectively. The two quantities are of course proportional, 

S = rSj^ (15) 

where the constant of proportionality r takes into account the exact shape of the UV 
background spectrum. If we express the UV background spectrum in terms of a broken 
power law, then from Eq. |^, we obtain 

= — ^r^- (16) 

c"Heii -J + Phi 

However, the latest models of the UV background spectrum depart significantly from a 
broken power-law. For example. Figure |l^ shows the UV background spectrum predicted 
by Madau, Haardt & Rees (MaHR; 1999). We have numerically evaluated the integral given 
by Eq. ||, using values for J(i/) kindly provided to us by F. Haardt. We obtain 5* ~ 230 
and S'l — 30, so that rMaHR — 8.8. [I'ardal et al. 199^ (see their Figure 6) present another 
model of the UV background spectrum resulting from their source model Ql with spectral 
index as = 1.8, stellar contribution, absorption model A2 and taking cloud re-emission into 
account. The stellar contribution is fixed at 1 Ryd to have an emissivity twice that of the 
quasars and is considerably softer than the MaHR's model or other models without stellar 
contribution. This Fardal et al. model has S ^ 930, Sl — 140 and rpcs — 6.9. 



We can set constraints on the UV background shape and intensity from the spectrum of 
Q 0302-003 for several reasons. First, the QSO is bright so that the proximity effect extends 
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Fig. 10. — Comparison between the Madau, Haardt and Rees (1999) spectrum of the UV 
background a.t z — 3 (left panel) and broken power laws of index /3hi = /^Heii — 1 (middle 
panel) and /3hi = 1, Pueii — 5 (right panel). For the middle and right panels, the values of 

JhcII and Jhi are inferred from the photo-ionization rates necessary to explain the proximity 
effect and the observed opacity in the He II opacity trough, respectively. The Madau et 
al. (1999) UV background spectrum has been multiplied by 1.47 to match our choice of 
cosmological parameters. 



- 31 - 



far away from it. Second, the largest A^hi cloud in front of it has a logA^ni = 15.4, so that 
even with values of 77 ~ 500, very few clouds have logA'neii > 17 where other radiative 
transfer effects start to be important. There is no 2;abs — ^em system nearly as strong in H I 
absorption as in HE 2347-4342 ( [Reimers et al. 1997 ), which leads to continuous absorption 



shortward of the He II edge for that object. The closest major absorber to Q 0302-003 at 
z = 3.2673 shows only a single C IV component, two O VI components and a H I column 
density of log A^hi = 14.66. There are no higher column density systems at 2; > 3.221. 

Third, the STIS spectrum covers a wide range in redshift, probing both the proximity 
region and regions far from the QSO. Without the proximity effect, the observations 
could not constrain the UV background shape and intensity independently. Rather, the 
observations yield the softness parameter, S, which is a function of both Jm/ Jueii, and the 
continuum slopes, This ambiguity is also shown in the middle and righthand panels 
of Figure [l^. In both panels, the UV background has a softness 5* = 800, but in the 
righthand panel, JhcII is two times higher which is compensated for by the steeper Lyman 
continuum slope. With the observation of the proximity effect, however, we gain a direct 
measurement of Tj^^jj by itself (not the ratio ^m/^Keii)^ which can then be combined with 
the observationally derived quantities to calculate Jneii and Jhi. 



5.3. Sensitivity of the models to input parameters 

Before presenting the results of our modeling, we first summarize the sensitivity of 
the model spectra to input assumptions or parameters, such as the presence/absence of a 
diffuse gas component ( §5.3.11 ), the UV background H I photo-ionization rate Fj^j, and its 



spectral shape as described by the softness parameter S. As shown in Eq. S is related 
to the conventional softness parameter, S'l = Jm/ Jneii, but we work in terms of S, because 
it directly takes into account the shape of the UV background spectrum. In § |5.3.2| , we will 



show that far from the QSO, the He II Lya spectrum depends only on S, but that close to 
the QSO, it depends on both F^j and S. Finally, we will also briefly discuss the effect of 
turbulent vs. Doppler broadening of the He II lines in § ^.3.3| . 

For simplicity, we assume that the slope of the QSO Lyman continuum evaluated 
near the He II Lyman-limit is the same as that of the H I Lyman continuum, i.e. 
'^ffl*^ = a§g^ = 1.9. We will then discuss the sensitivity of the results to a§gn in §|5.3.4 
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5.3.1. Presence of a diffuse gas component 

Diffuse gas in the form of numerous clouds of low column density (A^hi < 10^'^ cm~^) 
produces a decrease in the transmitted flux in the H I spectrum compared to the 'observed' 
line list, which is only complete down to A^hi < 10^'^ cm~^. This diffuse gas is especially 
important for the He II Lya forest at z ~ 3, where the optical depth is larger than unity 
for large values of the softness parameter. 



5.3.2. Dependence of rj on Fpjj and S 

Far from the QSO, where the second term of Eq. |12| is negligible compared to the first 
one, the value of the ratio of column densities r] = A^Heii/-^Hi is only proportional to S, 
through Eq. |^. This region of the spectrum can be used to constrain the spectral shape of 
the UV background. 

Close to the QSO, the second term in Eq. ^ becomes important. For reference, the 
H I photo-ionization rates due to the UV background and the QSO are equal ai z = 3.268 
if log Jhi = —21, = 1; and for the QSO H I Lyman limit luminosity shown in Fig. |. 
The same will be true for the He II photo-ionization rates if S'l = 4"*^^° = 13.9 (and if 
/^Hi — Pueii — Decreasing F^jj while holding constant the softness parameter, S, extends 
the redshift range over which the QSO has a significant effect on both H I and He II. 
Increasing 5* while holding Fj/jj constant has the same effect but only for He II, since this 
manipulation actually corresponds to decreasing F^^jj. Since decreasing Fjjj and increasing 
S have the same effect on the He II spectrum, the proximity effect can constrain the softness 
parameter if Ffjj is known. Similarly, limits on Fjjj can be derived from 5* which can be 
estimated from the run of r] far from the QSO. 



5.3.3. Doppler broadening 

As mentioned in § pl2| , we only consider turbulent line broadening (&Heii/&Hi = 1) to 
calculate the He II absorption spectrum. If instead, thermal motions were the dominant 
broadening mechanism (&Heii/^Hi = 0.5), then the He II opacity would be lower, and the 
UV background would have to be even softer in order to reproduce the observations. 
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5.3.4- Sensitivity of results to the QSO flux distribution 

The value of = 1.9 is secure, but stating that a§en =1-9 resuhs from an 

extrapolation. To our knowledge, the He II Lyman continuum region has never been 



observed, even in low-redshift quasars. As noted in § |5.1.1| , a softer QSO spectrum is 
theoretically possible but improbable since high-ionization absorption lines are detected 
in quasar spectra. If the spectrum were as soft as given by a§gj° = 11, then the QSO 
contribution to the photo-ionization rate of a cloud would be three times lower than for 
•^HeH — 1-9; and the proximity zone would be smaller. In order to match the observations, 
^Reii would have to be smaller, and the softness parameter 5* increased by a factor of 3. 



5.4. Results 

We now present the results of modeling He II Gunn-Peterson absorption in the 
spectrum of Q 0302-003. Comparing the modeled spectra with the observations allows us: 
(1) to confirm the absence of a diffuse gas component close to the quasar (§ |5.4.1j ); (2) to 



constrain the values of Jhi and S from the proximity effect (§ |5.4.2| ); (3) to show that diffuse 
gas component must be present far from the quasar, and to set limits on rj and S (§ |5.4.3| ); 
(4) to show evidence for a soft UV background ( $.4.4|) , which allows the UV background 



flux at the H I Lyman limit Jhi to be determined ( §5.4.51) . 

In practice, the calculations were performed using /S^i = /^Heii = 1 ^^e slopes of the 
UV background just blueward of the Lyman limits, and we present the models in terms of 
Jhi and S for easy comparison with other studies. Since the the fundamental variables are 
the photo-ionization rates, an easy conversion to Jj for other values of (5- can be done using 
Eq. |. 



5.4.1. No diffuse gas component close to Q 0302-003 

We explored a wide range in parameter space (Jhi, S) and found no way to reproduce 
the observed quasar spectrum in the proximity zone (1280 - 1300 A) if we use the combined 
line list., i.e. if we allow the presence of a diffuse gas component. Even for a UV background 
spectrum that is harder than the QSO's (cf. Figure 0), the 1297 A edge is too depressed, 
and the region between 1285 A and 1297 A is too fiat compared to the observed spectrum. 
On the other hand, we can achieve a satisfactory fit if we use only the observed line list. 
We therefore agree with Hogan et al. (1997) that there is no diffuse component close to the 
QSO. 
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Fig. 11. — Attempt to reproduce the observed spectrum (solid) using the combined 
Une Ust (dotted: Jhi = 1 x 10^^^ erg s~^ cm~^ Hz~^ sr~^, dashed: Jhi = 1 x 
10~^^ erg s~^ cm~^ Hz~^ sr~^). The 'best' match is obtained with S = 40, as shown here. 
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5.4-2. Constraint on Jhi and S from the proximity effect 

Here, we focus on reproducing the STIS spectrum in the proximity zone of the QSO, 
i.e. A = 1282 — 1301 A. A large number of observed pixel fluxes can be reproduced by a 
wide range of values of Jhi and S and thus are not useful for their precise determination. 
There are, however, two regions of ~ 5-pixel width, centered at ~1285 A and ~1290 A, 
that are quite sensitive to Jhi and S. Indeed, in the absence of a diffuse component, the 
predicted values of A^Heii for lines in these regions imply a rneii ~ 1 in the cores of the 
lines for the values of F^j and S of interest. Hence, these regions provide the best way 
to estimate Jhi and S. The opacity in the 1285 A region is somewhat larger than in the 
1290 A region, so that satisfactory fits are only obtained with a relatively small set of 
parameters. However, it is possible that the 1285 A region is affected by an ionizing source 
that could be also responsible for the Dobrzycki-Bechtold void seen in H I (see below). If 
so, the 1290 A region provides the stronger constraint on the UV background spectrum. 
Finally, we note that there are other spectral intervals within the proximity zone that have 
line-core opacities close to 1, but at the STIS resolution, these intervals are contaminated by 
strong lines nearby. Spectra at higher resolution are needed to better define these spectral 
intervals with r ~ 1 and thus allow a better determination of F^j and S. 



Figure ^ compares two model spectra with the normalized observed spectrum. The 
two models, with their different values of S, bracket a best model of the proximity effect. 
If Jhi were fixed at a lower value, e.g. 1 x 10~^^ erg s~^ cm~^ Hz~^ sr^^, then a very large 
value of S* 4000 would be required. On the other hand, the background radiation field 
could be as high as Jm = 1 x 10"^^ erg s"^ cm"^ Hz"^ sr"^ if 400 < S < 800. 



5.4.3. Determination of S far from the quasar 

We now estimate the softness parameter S at lower redshifts. It quickly appears that a 
diffuse component is needed to reproduce the high He II opacity in the D-B region. We thus 
created a line list that contains only the observed lines at z > 3.220 (i.e. in the proximity 
zone) and all the lines from the 'combined' line list at z < 3.220. This boundary redshift 
corresponds to the redshift of the strongest H I absorber in the spectrum of Q 0302-003. As 
this absorber is relatively close to the quasar, it efficiently blocks He II ionizing flux from 
the quasar. 



Figure |T3| compares the observed spectrum with model spectra computed for two 
different values of the softness parameter S. It shows that a hard UV background with 
S ~ 120 produces a good match to the He II opacity gap at A = 1230 A {z = 3.05), and it 
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Fig. 12. — Models of the proximity effect region for two values of the softness parameter 
^ — rRi/rHeii (i^ gr^y)- The sohd line represents the normalized STIS spectrum. 
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provides a reasonable fit to the low-redshift regions of tlie spectrum, 1150 A < A < 1175 A. 
Tlie region, 1205 A < A < 1210 A, needs a somewliat softer background witli S ~ 300. 
However, tfie UV background must be quite soft, i.e. 5* = 800 — 1000, to explain the 
observed opacity over most of the STIS spectrum including the Dobrzycki-Bechtold 'void' 
(the D-B region). An even softer UV background produces too large an opacity: if 
S ~ 1400, the model produces only 1/6 of the flux observed in Region A. 

We stress that the method here automatically takes into account changes in the gas 
density traced by the Lya clouds, as long as most of the hydrogen is ionized and helium is 
doubly ionized. Consequently, the features that our model spectra cannot reproduce with a 
unique value of the softness parameter are likely due to local changes in the ionizing flux. 
For example, the fact that the strong He II opacity in the Dobrzycki-Bechtold 'void' can 
be reproduced by a diffuse component suggests a nearby source that is able to completely 
ionize the hydrogen but unable to doubly ionize the helium. Such a source must have a 
much softer spectrum than a typical QSO. Other examples include the low He II opacity 
at 1230 A {z = 3.05) which corresponds to a local decrease in the Jui/Jucu ratio, and 
suggests an AGN close by, or some overlapping Stromgren spheres of more distant sources. 
The marginal decrease in the He II opacity at 1207 A may be similarly explained. In 
addition, some regions of low He II opacity like the one at ~ 1182 A are expected in our 
model; we discuss them in § |6.2[ Finally, the extent of the region 1150 A < A < 1175 A and 
the observation of HS 1700+6416 (Davidsen et al. 1996) suggest that the UV background 
presents a harder spectrum at 2; < 3 at which the Universe is effectively transparent. 

Figure [1^ shows how the ratio of column densities t] = A^Heii/A^Hi varies as a function 
of wavelength (redshift). All the curves assume the same intensity of the UV background, 
log Jhi = —20.60, and use the same line list (observed line list in the proximity zone, 
combined line list further away from the QSO). They differ only in the softness of the UV 
background. In all cases, 77 increases gradually with distance from the QSO. A similar figure 
for larger Jhi would display flatter curves. In all the cases which reproduce the data well, 
rj > 350 far from the quasar, as shown by the top curve. 

5.4-4- Evidence for a soft UV background 

From the previous sections, and as also inferred from previous observations of the 
He II Gunn-Peterson effect, we find little support for a 'hard' (Sl ~ 30) UV background 
as predicted by Haardt & Madau (1996) and Madau, Haardt & Rees (1999). Instead, in 
addition to the direct measurement of theii described in section H, our modeling shows 
evidence for a soft background both from the proximity effect and from the redshift domain 
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Fig. 13. — Model vs observed spectrum. The thick sohd hne represents the observed 
spectrum. The gray regions shows the model spectra: above: with S = 120 that best 
matches the low opacity region at ^ = 3.05 and produces an acceptable fit in the regions 
1150 A < A < 1175 A; below: with S = 800 (and log Jhi = —20.6) giving a good fit over 
most of the observed spectral range. 
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Fig. 14. — Behaviour of rj as derived by our model using the 'observed' hue hst in the 
proximity region and the 'combined' hne hst far from the quasar. Three different values 
of S are presented for log Jhi = —20.60, although only the S = 800 case models the data 
well. Note that lower values for Jhi are allowed, but require even larger values for S, so that 
rj ^ 200 far from the quasar in all the cases which reproduce the data acceptably. 
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away from the quasar ionizing region. The value of the softness parameter S ~ 800 is close 
to the value ~ 930 derived for the Fardal et al. (1998) model (their Figure 6, source model 
Ql with spectral index as = 1.8, stellar contribution fixed at 1 Ryd to have an emissivity 
twice that of the quasars, absorption model A2 and taking cloud re-emission into account), 
which suggests that a significant stellar contribution is required. 

In order to compare the ratio = Jm/ JhcU of H I to He II UV background intensities 
with other values found in the literature, we have to assume a shape for the UV background 
spectrum. For a broken power-law model with /?hi = Pneii = 1, the proportionality constant 
r = 4 (Eq. |I6]) so that S'l = 200. Instead, if we use the shape of the spectrum obtained by 



Fardal et al. (1998), r = rpcs - 6.9 (§§g) and 5l ~ 120. 



5.4-5. A determination of Jm 

Since a softness parameter S ^ 800 is required by our models over most of the range of 
the STIS spectrum, the constraint imposed by the proximity effect sets log Jhi = —20.6 for 
/3m = 1 (cf. Fig. m. 



Taking into account the differences in the assumed cosmological parameters and UV 
background spectral slopes between the ones we assume and those chosen in the following 
studies, this value of the UV background is smaller than the 3a upper limit derived by a 
recent attempt to detect Lya emission in Lyman-limit absorption systems ( [Bunker et al.| 
1998 ), but is a factor of 2 larger than the one obtained by [Scott et al. (199^ and a factor of 
~ 3 larger than the one derived by [Giallongo et al. (1996)| from the study of the Lya forest. 



6. THE OPACITY GAP AT z = 3.05 

We now return to the windows of transmitted flux seen in all three sightlines probed 
by STIS (cf. Figure These windows give insight into the reionization history of the IGM 
and how the reionization of He occurred. Indeed, their observed properties provide a direct 
point of comparison with the growing number of theoretical studies that predict the He II 
opacity with respect to the UV background flux, IGM density and other physical aspects of 
the early universe ([Croft et al. 1997| ; [Fardal et al. 199"8| ; [Zheng fc Davidsen 1995[ ; [Zhang et 
al. 1998[ ; MHR and references therein). 



The results of our models described in §5 imply that the opacity gap at z = 3.05 is 
most likely produced by a local source able to doubly ionize helium, or the result of a 
change in the Jm/ Jneii ratio (due, for example, to additional ionizing sources or to increased 
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transparency of the medium to He II Lyman continuum radiation). In this section, we 
analyze this opacity gap in more detail. We first describe in §6.1 the observed properties 
of these regions of enhanced transmitted flux seen in the Q 0302-003 spectrum. We then 
examine three proposed explanations of He II opacity gaps. In §3^, we consider the 
possibility that low opacities are the signatures of voids, or underdense regions [p/p ~ 10~^). 
For our second interpretation, we explore in §^]3| the hypothesis that opacity gaps are 
associated with regions in which helium has been doubly ionized by nearby discrete sources. 
Here, we have used the estimated UV background from §^ to test whether an AGN is 
required near the z = 3.05 gap. Finally, in §6.4, we briefly examine the possibility that 
opacity gaps are regions that have been collisionally ionized by shock-heated gas. 



6.1. Observed properties of He II opacity gaps 

Table ^ gives some of the measured characteristics of the observed opacity gaps in all 
three sightlines probed by HST. Successive columns list the redshift of the opacity gap and 
its comoving distance from the QSO (AD), the size of the gap (in A, km s~^, and comoving 
Mpc), and finally, the HST instrument used to obtain the observation. For consistency, 
the redshift and widths of the two gaps in the spectrum of HE 2347-4342 were measured 
from GHRS data rebinned to the (lower) resolution of STIS. The gap size is defined as 
the interval of continuous pixels with flux > la, except for PKS 1935-692, where the flux 
redward of the gap is not zero. In that case, the boundary of the gap was taken as the 
wavelength at which the flux equals the mean flux over 1248-1264 A. 

Since all three lines of sight observed by HST show opacity gaps, we assume that the 
line of sight to Q 0302-003 is typical. We can therefore calculate the filling factor of the 
region responsible for the z = 3.05 opacity gap as the ratio of the path-length through 
the gap, Ai ~ 17.2 comoving Mpc, with respect to the total path-length probed by the 
STIS spectrum, AL. We obtain AL = 395 comoving Mpc if we integrate over the observed 
Gunn-Peterson trough (2.7609 < z < 2.9879 and 3.0142 < z < 3.2183; that is, avoiding the 
proximity-effect region and the region contaminated by geo-coronal Lya), or AL = 180 
comoving Mpc, if we only integrate down to z = 3.0142, just longward of geo-coronal Lya. 
The resulting filling factors are 0.04 and 0.10 respectively. If we include all three sightlines, 
which probe a total of AL ~ 560 comoving Mpc, then we obtain a filling factor of 0.09. If 
we assume that each region of enhanced transmission is a sphere whose diameter is given 
by the average size of the opacity gaps, d = 13 comoving Mpc, then the volume probed by 
the spectra is vr {d/2Y AL ~ 7.4 x 10^ comoving Mpc^, and the number density of opacity 
gaps is ngap ~ 5 x 10~^ Mpc"^. 



Table 3. Hell Absorption Gaps 



Sightline Zgap AD Linear extent Instrument 

(comoving Mpc) (A) (km s"-*-) (comoving Mpc) 



Q 0302-003 3.052 

PKS 1935-692 3.100 

HE 2347-4342 2.865 

HE 2347-4342 2.814 



202.5 5.8 1420 

69.8 4.7 1130 
19.0 2.3 590 

67.9 4.0 1043 



17.2 STIS 

13.6 STIS 

7.2 GHRS 

12.8 GHRS 



Table 4. H I, metal absorbers near the z = 3.0526 He II opacity gap 



ion 



b (km s ''") 



log A^ion {cm-y 



CIV 3.046980 ± 0.000018 6.2 ±2.2 12.22 ±0.11 

CIV 3.047275 ± 0.000015 11.4 ± 2.0 12.62 ± 0.05 

HI 3.045994 ± 0.000039 28.1 ±2.6 14.50 ± 0.09 

HI 3.047073 ± 0.000036 63.8 ± 2.0 15.32 ± 0.02 

HI 3.048479 ± 0.000039 35.6 ± 2.0 14.29 ± 0.06 



^Upper limits (4ct, assuming z = 3.0472) are: C III A977 < 11.45, C II A1335 < 11.93, 
Si IV A1393 < 11.35, Si III A1206 < 11.05, Si II A1526 < 12.06, N V A1238 < 11.93, 
O VI A1031 < 12.26. 
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6.2. He II opacity gaps as low-density regions 

Cosmological hydro-dynamical simulations predict growing density fluctuations with 
time, and MHR interpret opacity gaps in the He II Gunn-Peterson absorption trough as the 
spectral signatures of underdense regions. Below, we test their interpretation by comparing 
the predicted properties of low-density regions to the observed widths and amplitudes of 
opacity gaps. 

Sizes of underdense regions. MHR modeled the effects of a clumpy IGM on the 
reionization history of helium. They concluded that there were sufficient photons to ionize 
most of the helium by z = 3 unless very luminous QSOs were the only sources of He II 
ionizing photons. Based on the He II Lya opacities we measure, they estimate that at 
z = 3, helium is doubly ionized up to an overdensity of p/p ~ 12, with a volume fraction 
Fy ~ 0.9943 of helium doubly ionized, and a mean free path of 4-Ryd photons of ~ 1500 
km s~^. Fluctuations in the density and ionizing background should produce opacity 
fluctuations with a similar scalelength. The predicted mean free path is consistent with the 
largest of the He II opacity gaps (cf. Table 3). 



Transmission of underdense regions. In § |5.4.3| , we noted that gaps in the He II opacity 



appear in our simulations due to random fluctuations in the diffuse IGM (i.e. logiVni < 13). 
An example of such an opacity gap can be seen at A ~ 1182 K {z = 2.890) in Fig. 0. 
Could similar fluctuations produce the prominent z = 3.05 opacity gap, which has optical 
depths as low as r = 0.5? To answer this question, we obtained 100 realizations of the 
simulations described in §|| using S = 800 and log Jhi = —20.6. As mentioned in § ^.4.2| , this 



value of S can account for the large opacities in the Gunn-Peterson trough; lower values of 
Jhi lead to larger 5* values. We then measured the distribution of opacity gaps produced 
by our simulations. We found that none of our 100 simulations produced a gap having 
an optical depth smaller than r = 1.3 over the whole spectral range covered by our STIS 
spectrum. Instead, the simulations produced an average of 1.2 opacity gaps with optical 
depths r < 2.0. Most occur at low redshift where the number density H I Lya lines is low, 
or less frequently, near the redshift of the QSO where excess ionization by the QSO is not 
negligible. In contrast, the observed opacity gap at z = 3.05 is located in a region unlikely 
to show an opacity gap caused by fluctuations in the diffuse IGM. In fact, only one of our 
simulations produced a gap having an optical depth smaller than r = 1.8 in the redshift 
range 3.0 < z < 3.1. 

In summary, the predicted sizes of underdense regions are compatible with the observed 
widths of the opacity gaps. However, our simulations indicate that the underdense regions 
cannot produce the distribution and amplitude of the observed transmission. This failure 
constitutes a strong argument against the hypothesis that the z = 3.05 gap arises as a 
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fluctuation in the diffuse gas component of the IGM. 



6.3. He II opacity gaps produced by discrete ionizing sources 

Another possibihty advanced by Reimers et al. (1997) is that He II opacity gaps may 
be caused by discrete ionizing sources along or near the QSO hne of sight. In this case, 
the UV source would most likely be an AGN or QSO since star-forming galaxies will not 
substantially ionize intergalactic He II ( [bardal et al. 1998 ). Support for a discrete ionizing 



source comes from the optical Keck HIRES data of Q 0302-003. As shown in Figure |T3|, 
the Keck spectrum reveals eight C IV absorption complexes in the redshift range probed by 
STIS. The redshifts of these complexes are indicated by vertical bars. We follow Steidel 
(1990) in defining a complex as a group of systems spanning less than 1000 km s~^ in 
velocity space; in each of the relevant complexes towards Q 0302-003, the components span 
less than 200 km s~^. Two of the complexes, at z = 3.22 and 3.27, are at the edge of (or 
within) the proximity zone of the QSO. One C IV complex, at z = 3.0, falls too close to 
the H I Lja geo-coronal emission line to judge whether it has associated transmission of 
He II Lya. Of the remaining five, four C IV absorption complexes, at z = 2.79, 2.83, 2.97 
and 3.05, are close to regions with diminished He II Lya opacity. Assuming a Poissonian 
distribution of C IV systems, we estimate a probability, p = 12% of one coincidence of a 
C IV absorber with a He II opacity gap; four coincidences are much less likely {p = 0.02%). 

The window at z = 3.05 is of course the most prominent of the three. Figure |16| 
shows detailed plots of the C IV and H I absorption by this system. Table 4 gives the 
detailed measurements of absorbing systems near the major He II opacity gap at z = 3.05. 
Successive columns list the ion (and spectral line), measured redshift, Doppler parameter b, 
and column density. The table also gives Aa detection limits for a series of C, Si, N and O 
ions. 

Below, we test the hypothesis that the C IV absorber at z = 3.05 could be photo-ionized 
by a hypothetical AGN that is also responsible for the adjacent He II opacity gap. Our 
diagnostics include: the required luminosity of the ionizing source, its spectral energy 
distribution and the space density of galaxies or AGN at 2; ~ 3. We shall find consistent 
support for the discrete ionizing source hypothesis. 

Luminosity requirements. To first approximation, the luminosity of the putative 
He"'"-ionizing source relative to the QSO scales with the relative volume of their Stromgren 
spheres. The radius of the QSO Stromgren sphere is given by the physical extent (luminosity 
distance in the local frame) of the proximity effect, A^prox ~ 14.4 Mpc, while the radius of 
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Fig. 15.— Location of CIV absorbers toward Q 0302-003. The STIS spectrum of Q 0302-003, 
binned by 3 pixels for clarity, is shown by the thick solid line. The optical HIRES spectrum 
(thin solid line) has been divided by a continuum fit to remove the Lja emission line, scaled 
by 2.5 X 10^^^ to match the STIS spectrum redward of 1300 A, and offset by 1 x 10~^^ for 
clarity. The HIRES wavelengths have been divided by the ratio of the wavelengths of Hell 
to HI hya. The vertical ticks show the redshifts of detected C IV absorption systems. Flux 
is in erg cm~^ s""^ A~^. 
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Fig. 16. — The Voigt profile fit in velocity space, in units of normafized fiux vs. wavelength 
in A, to the z = 3.047 C IV doublet toward Q 0302-003. The H I Lya and Ly/3 regions 
required a three component fit, while the C IV doublet required only two components. 
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the ionized region caused by the putative ionizing source is the half-width of the opacity 
gap A£gap ~ 2.15 Mpc. Given these size estimates, the hypothetical source would need to 
be only 0.3% as luminous as the QSO (V = 17.4, |Veron-Cetty fc Veron, 195BD and would 



have an apparent magnitude of V = 23.4 if it were on the line-of-sight to the QSO and 
radiates isotropically. A more detailed luminosity estimate with Kraemer's (1985, |Kraemer| 
et al. 19"9^ ) photo-ionization code supports the conclusion that a bright Seyfert I galaxy 



about 0.3% as luminous as Q 0302-003 could produce the observed gap at z = 3.05 if it 
were directly along the line of sight to the QSO. A source offset from the line of sight would 
need to be more luminous. 

The ionizing flux distribution near the C IV absorber. Let us assume that the metal-line 
system aX z = 3.047 originates in the halo of the host galaxy (or a nearby galaxy in the host 
cluster) of an AGN, the putative source of the ionization which produces the opacity gap. 
To explore whether this C IV absorber could be photo-ionized by the hypothetical AGN, 
we constructed photo-ionization models with CLOUDY ([Ferland et al. 1998|) . For the 



model calculation, we assumed the usual plane-parallel geometry and allowed the gas to be 
ionized by a model AGN spectral energy distribution (SED) and/or a Fardal et al. (1998) 
UV background at z = 3.0. The Fardal et al. UV background spectra are more consistent 
with the results derived in than the Haardt & Madau (1996) or Madau, Haardt & Rees 
(1999) ones because of their larger values for = Jm/ Jneii- We tried the Fardal et al. 
backgrounds due to QSOs with and without a contribution from stars (see their Figure 
6). We used Ferland et al.'s (1996) model of an AGN SED, where the "big blue bump" is 
approximated as a power law with a UV exponential cutoff with characteristic temperature 
Tbb- We calculated the C IV/Si III, C IV/Si IV, C IV/N V, and C IV/0 VI column density 
ratios as a function of the ionization parameter, U = ra^/nn, where is the H I ionizing 
photon density and wh is the total hydrogen number density. 

Figure |1^ shows a sample CLOUDY calculation for the case in which the ionizing flux 
is dominated by the UV background, including a contribution from stars. We find that the 
model is in agreement with the observed constraints (Table 4) for 10~^'^° U < 10~°'^^, 
which corresponds to 1.5x10"^ < n-n/ J21 < 5.3x10"^ cm~^ where J21 = Jm/lO"^^ erg 
s~^ cm~^ Hz~^ sr^^. Since these are reasonable particle densities for the distant outer halo 
of a galaxy, we conclude that it is not necessary to have an AGN nearby to explain the 
observed high level of ionization. However, we find from another CLOUDY calculation 
in which the ionizing flux is due to a nearby AGN that the model column density ratios 
can be brought into agreement with the observations over a reasonable range of U as long 
as the characteristic temperature of the "big bump" Tbb is less than ~ 10^'^ K. Since 
luminous QSOs tend to have higher values of Tbb (Hamann 1997), this indicates that the 
putative AGN must be at least a low-luminosity QSO. We conclude that it is possible, but 
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Fig. 17. — Column densities of high-ionization species, predicted by a photo- ionized gas 
model for the C IV absorber at Zabs = 3.047 as a function of the ionization parameter U 
(bottom axis) and total hydrogen number density nn (top axis). This model assumes that 
the ionization is predominantly due to the UV background from QSOs and AGN and adopts 
the radiation field calculated by Fardal et al. (1998) for z = 3.0. The relative heavy element 
abundances are assumed to be solar, and the overall metallicity [M/H] = -2.5 with log Nm 
= 15.32. 
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not necessary, that a low-luminosity AGN near the sight line at z ^ 3.05 ionizes the region 
around the He II opacity gap and may also play an important role in the ionization of the 
nearby C IV system. 

The space density of opacity gaps vs. galaxies. How plausible is it that regions 
of low He II opacity are produced by galaxies, or more likely, galaxies associated with 
AGN? To answer this question, we compare the number density of He II opacity gaps, 
'^gap ~ 5 X 10~^ Mpc~^, with the number density of star- forming galaxies and AGN as 
revealed in narrow-band or emission-line surveys of galaxies at high-redshift. These surveys, 
which are summarized in Table 5, show that the comoving space density of galaxies at 
2.4 < z < 3.4 is ngai ~ 60 — 410 x 10~^ Mpc~^. Hence, the space density of opacity gaps 
is about 1-8 % that of galaxies. If the AGN fraction among galaxies were similar, it would 
be consistent with the hypothesis that He II gaps are produced by such discrete ionizing 
sources. The incidence of low- luminosity AGN at high redshift is poorly known. About 
10% of Lyman limit dropout galaxies at z ~ 3 are QSOs or AGN (|Steidel et al., 19961) with 
a space density of ~ 1 x 10~^ comoving Mpc~^. An independent estimate by Teplitz et 
al. (1998) from an extrapolation of the QSO luminosity function gives a similar number 
density at 2.3 < z < 2.5 for an Ha emission line flux of 10~^^ erg cm~^ s^^. This value 
is within a factor of two of the number density of He II low-opacity regions, so it is quite 
plausible that the gaps are produced by AGN. 

6.4. Can the gaps be caused by shock-heated gas? 

In hydro-dynamical cosmological simulations, galaxies are found in the vicinity of hot, 
collisionally ionized gas regions. This is not only because the collapse of initial density 
perturbations leads to shock heating of the gas; it also enables the condensation of cool 
objects to form galaxies (Dave et al. 1999; Gen Sz Ostriker 1999). It is likely that the C IV 
absorber near the He II opacity gap at z = 3.05 arises in a galaxy halo. Thus, it is possible 
that the low He II opacity is basically due to the high temperature of the region in which 
the galaxy is located. The volume fraction of gas with temperature 10^ < T < 10'' K is 
roughly 0.01 at z ~ 3 (see Figure 2 in Gen & Ostriker 1999), which is a factor of 4-10 times 
lower than that of the observed opacity gaps. It may be possible to explain the discrepancy 
by uncertainties in the hot gas fraction as computed by the hydro-dynamical codes; the 
matter clearly deserves further study. 
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6.5. Summary 

In summary, we have considered three explanations of the high-transmission regions 
in the He II Gunn-Peterson absorption trough: (a) they result from underdense regions 
in a clumpy IGM; (b) they are regions ionized by nearby AGN; and (c) they arise in hot, 
collisionally ionized regions. We find that the observed opacity gap at 2; = 3.05 is located in 
a region unlikely to show an opacity gap caused by fluctuations in the diffuse IGM. Discrete 
ionizing sources, however, are consistent with the observed gaps for the following reasons: 
(1) The luminosity required is consistent with an AGN. (2) There is a C IV absorber at 
a velocity separation of Aw = 350 km s~^. (3) The ionization state of the C IV absorber 
does not preclude the existence of a nearby AGN, although it does rule out the possibility 
of a bright QSO in the immediate vicinity. (4) The space density of He II opacity gaps is 
consistent with the best available estimates of AGN at z ~ 3. Collisionally ionized regions, 
which are predicted in galaxy formation models, do not produce a high enough filling factor 
of ionized gas in current models. 



7. DISCUSSION AND CONCLUSIONS 

Our analysis of the STIS spectrum of Q 0302-003 leads to the following conclusions. 

1. Opacity of He II Lya at z ^ 3. Away from distinct opacity gaps that we could identify 
in the He II absorption, we made direct measurements of the average transmission of the 
IGM. We found values that ranged from / = 0.0086 for an interval centered at z = 3.15 
to / = 0.15 for another zone centered aX z = 2.82. If the opacity were constant within 
one of these intervals (which it is not), these intensities correspond to r = 4.75 and 1.88, 
respectively. A stronger statement can be made about the IGM opacity when we note that 
there are no regions near 2; = 3.15 that have a relative transmission greater than 0.03. If 
we apply this result to the expectations for opacity fiuctuations computed by Fardal et 
al. (1998), we find that the average transmission must be even lower than what we could 
measure, with representative values of optical depth greater than 6.0, depending on some 
initial assumptions. 

2. Evolution of He II T(Lya). The total He II opacity rises more rapidly with redshift 
than previously thought (e.g. [Fardal et al. 1998| , Zhang et al. 1998). The observed data 



are compatible with an opacity break occurring between z = 2.9 and z = 3.0 as suggested 
by other lines of evidence presented by Songaila (1998). This rapid rise in opacity, however, 
is overshadowed by the fact that at 2; ~ 3, He II is such a minor fraction of helium 
(~ 1.6 X 10~^). Clearly, reionization of helium is essentially complete by 2; ~ 3. 
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3. Comparison with cosmological models. Observations of He II Lyct absorption give a 
direct connection with cosmological simulations if differences in resolution are taken into 
account. A comparison of Zhang et al.'s model to the observations indicates that their 
predicted He II Lya opacity is 3-5 times lower than observed by STIS. 

4- The UV background and the He II Lya opacity. We constructed detailed models of the 
Gunn- Peterson trough over the observed range in redshift. In agreement with Hogan et al. 
(1997), we find that the observed proximity effect cannot be reproduced in the presence 
of a diffuse component. Instead, a model based only on known H I Lya lines provides a 
good match to the observations. The model also requires that the UV background has a 
soft spectrum with a softness parameter S = T'^i/T^^n certainly larger than 400 and close 
to 800, indicating a significant stellar contribution. We can also set the constraint that if 
Jm ^ 10~^^ erg s'^ cm''^ Hz"^ sr^^ then S > 4000. 

Outside of the proximity region, most of the He II Gunn-Peterson trough requires a 
large value of 77 = -/VHeii/-/VHi — 350 ±50, which in turn requires both a diffuse gas component 
and a soft UV background with values of 5" ~ 800. In particular, these constraints are 
needed to explain the high He II opacity in the Dobrzycki-Bechtold void. As also deduced 
from the observations of HS 1700+6416 (Davidsen et al. 1996), our model indicates that 
the UV background has a harder spectrum at 2; < 3 as the Universe becomes transparent 
to He II Lyman continuum radiation. 

5. Nature of the opacity gaps. We considered three explanations for the He ll opacity 
gaps observed in QSO spectra: shock-heated gas, low density regions demarcated by gaps 
in the H I Lya forest, and regions ionized by a discrete, local sources. The shock-heated 
gas scenario, as calculated by current hydro-dynamical codes, is inconsistent with the 
observed opacity gap filling factor. The clumpy IGM model of Miralda-Escudc, Haehnelt 
& Rees (1999) predicts low-opacity regions with a mean free path similar to that of the 
observed He ll opacity gaps. We can reproduce the observed He ll opacity gaps only if 
77 ^ 100, contrary to other regions of the spectrum where rj ~ 350. The softness parameter 
corresponding to 77 ^ 100, S ^ 200, is consistent with a local hard ionizing source, such as 
an AGN. 

In the case of the z — 3.05 opacity gap, there are further supporting reasons for the 
presence of such an ionizing source. (1) The luminosity required is consistent with a bright 
Seyfert I galaxy. (2) There is a nearby C IV absorber at a velocity separation of Av — 350 
km s~^, possibly indicating the existence of a nearby AGN. Three other regions of lower 

He II opacity can also be associated with a C IV system. We estimate at p = 0.02% the 
probability that four opacity gaps are associated with C IV absorbers by chance. (3) The 
ionization state of the C IV absorber does not preclude the existence of a nearby AGN. (4) 
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The space density of He II gaps (5 x 10 ^ comoving Mpc ^) is consistent with the best 
available estimates of AGN at 2; ~ 3. 
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Table 5. High-Redshift Emission Line Galaxy Surveys 



redsiiitt 


density 


3(7 tiux 


iNo. 


Keterence 




(10-^ Mpc-=^)^ 


(10-^6 erg cm-2 s"^)'' 


found 




2.3-2.4, 0.89 


6 


4.8 


18 


Mannucci et al. 1998 


2.3-2.5 


41 


1.0 


5 


Teplitz et al. 1998 


3.4 


17 


0.09 


12 


Cowie & Hu 1998 



'^Comoving space density and star formation rate were calculated for our cosmological 

parameters. 

^Only objects with line fluxes greater than 1 x 10~^^ erg cm~^ s~^ were included in the 
tally. 
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